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Neutrinos at energies above TeV can serve as probes of the stellar progenitor and jet dynamics of 
gamma ray bursts arising from stellar core collapses. They can also probe collapses which do not lead 
to gamma-rays, which may be much more numerous. We calculate detailed neutrino spectra from 
shock accelerated protons in jets just below the outer stellar envelope, before their emergence. We 
present neutrino flux estimates from such pre-burst jets for two different massive stellar progenitor 
models. These should be distinguishable by IceCube, and we discuss the implications. 
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I. INTRODUCTION 

The gamma-ray bursts (GRB) which have so far been 
accurately localized are associated with regions of active 
star formation, and their progenitors are thought to be 
massive stars. The leading model for such bursts involves 
a relativistic jet, produced following the collapse of the 
core of the massive stellar progenitor 1]. In this model 
the 7-rays are produced by synchrotron or inverse Comp- 
ton radiation from Fermi accelerated electrons in opti- 
cally thin shocks (see for a review), after the jet has 
emerged from the stellar envelope. The same optically 
thin shocks should accelerate relativistic protons Q , and 
lead to ^ 100 TeV neutrinos via interactions with the 
observed MeV 7-rays ,4]. However, while the jets are 
still inside the star, shock-accelerated protons can pro- 
duce ^ TeV neutrinos through photomeson interactions 
with thermal X-rays in the sub-stellar jet cavity |^. 

In this paper we discuss a more general class of massive 
stellar collapses, in which jet formation may be ubiqui- 
tous, but not all of which emerge to be associated with 
detectable GRBs. Before their successful or failed emer- 
gence from the star, the jets can accelerate protons which 
undergo a more complex sequence of high energy inter- 
actions than previously realized. These depend not only 
on the jet and central engine characteristics but also on 
the location of the shocks and on the outer dimensions 
of the stellar progenitor, thus providing potentially use- 
ful diagnostics for the type of progenitor as well as the 
jet and shock parameters. Protons accelerated in sub- 
stellar jet shocks first undergo photomeson interactions 
with thermalized shock photons, as well as pp, pn interac- 
tions with thermal nuclcons in the jet frame. This mod- 
ifies the relativistic proton spectrum reaching the end of 
the jet cavity, where the protons undergo a second set 
of photomeson interactions with stellar X-ray photons 
and pp, pn interactions with cold nucleons in the stel- 
lar frame. The fraction of collapses producing jets which 
subsequently emerge from the star to produce electro- 
magnetically detectable GRBs are expected to be pre- 
ceded by a precursor neutrino signal at energies > TeV, 
which is significantly different from the previously cal- 
culated ^100 TeV neutrino signals coincident with the 



7-rays The fraction of stellar collapses leading to jets 
which do not emerge would have similar neutrino signals, 
but they could be more numerous and hence their diffuse 
flux could be more important. 

We discuss our jet models in Sec. El proton and elec- 
tron acceleration in the internal shocks in Sec. IIIII and 
proton interactions in Sec. IIVI We discuss neutrino pro- 
duction mechanisms in Sec. and calculate observed 
neutrino flux in Sec. IVII We summarize and discuss 
implications of our results in Sec. IVIII 



II. COLLAPSAR AND JET MODELS 

We take a simplifled model for the jet and progenitor 
star, similar to that described in Ref. The GRB 

progenitor is taken to be a massive star with a He core 
and H envelope. As an example, the parameters chosen 
are a core radius rjjo ~ 10^^'^ cm , where the density is 
Phc 10~'^p_3 g/cu^. We take two cases, one case (H) 
where the core is surrounded by an H envelope of size 
~ 10^^ cm, where the density pn ~ 10^''p-7 g/cm^, 
and another case (He) where a similar He core has lost 
its surrounding H envelope. Numerical simulations of 
coUapsar models leading to black hole formation (e.g. 
indicate that a relativistic jet can be launched along the 
stellar rotation axis, powered either by thermal uD an- 
nihilation or MHD stresses coupling to the black hole 
and to the debris disk falling back onto the hole). The 
jet life time is limited by the gas fall-back time onto the 
black hole, and should be comparable to GRB durations, 
< 10^ s. Variability in the efective jet Lorentz factor 
arises both from variability in the black hole accretion 
rate, and from instabilities at the jet-stellar interface on 
the way out from the star. 

Internal shocks occur at radii Tsh ~ c6tT^ in the highly 
variable jet outflow, where we expect fluctuations in the 
bulk Lorentz factor P 10^ — 10"^ over a wide range of 
time scales, 0.1ms 5t < \ s. The observed 7-rays are 
produced in jet internal shocks between regions with high 
enough P, 5t, which occur in the optically thin environ- 
ment outside the progenitor star (e.g. 0]). Fluctuations 
with lower P, 5t are also likely to produce internal shocks 
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at smaller radii below the stellar surface |^ , which will be 
optically thick. In observed 7-ray light curves, variability 
on sub-ms timescales is present, the power at high fre- 
quencies being subject to Poisson noise uncertainties 
due to low photon counts. In sub-stellar shocks with 
smaller lengthscales and Lorentz factors the high fre- 
quency power could be relatively more important. Such 
episodes of lower T or shorter variability could arise from 
instabilities leading to the up and down variation of F 
during propagation of the jet inside the star, as shown in 
numerical simulations 0. 

Here we are interested mainly in these opaque, sub- 
surface internal shocks. We consider the jets at a time 
when they are still inside the star, but have advanced to 
within a factor 1/few of the outer stellar surface. As spe- 
cific numerical examples we take the internal shock radii 
at VshM ^ al0^^-^ri2.5 cm and rsh,}ic ~ alO^rn cm, 
where a < 0.5, and the jet termination radii are taken to 
be rjet,H ~ 10^^-^ri2.5 cm and r^et, Ho lO^rn cm for 
the two types of progenitors. The internal shocks at Vsh 
provide a particle accelerator sufficiently far inside the jet 
cavity, where they are not subject to the large radiative 
losses in the high photon and nucleon density which are 
present further out near the jet termination radius rjet, 
i.e., we consider rsh < rjet < t^. The jets are optically 
thick to Thomson scattering and the density of photons 
from the jet termination surface drops off very fast with 
distance measured inwards into the jet . The shocks in 
the sub-stellar shocks are likely to be collisionless, since 
the typical collisionless width c/ujp^i (where tOp^i is the 
ion plasma frequency) is much smaller than the radiation 
scattering and the particle collision mean free paths. We 
introduce a parameter £op that characterizes the fraction 
of jet energy that is dissipated in the sub-stellar shocks 
(this energy is given back to the jet, due to the high opti- 
cal depth, except for the fraction that escapes as neutrino 
energy). The neutrino fluxes we derive will be propor- 
tional to £op- We cannot estimate Cop directly, since we 
do not see these collisions. However, it is likely to be 
significant due to the fluctuations in jet parameters, and 
we will assume Sop 1/2. 

For observed isotropic GRB luminosities Ll^° — 
10^^L52 ergs/s at the termination radius of a sub-surface 
jet, the jet-frame proton number density is 
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Equating the pressures behind the forward and reverse 
shocks one finds that the jet head moves at mildly rela- 
tivistic or subrelativistic speeds with Lorentz factor Q 



The termination shock where the jet head impacts the 
star heats up the stellar plasma. The temperature of 
thermalized shocked photons in the jet head frame, from 
^so ^ 4j,r%tcrlaBT\ is 



l-5Lli'rlli'r-\/!keY (H) 
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The corresponding number density of these photons in 
the jet head frame, which is approximately the same in 
the stellar plasma and in the observer's frame, is 

^ / 4.3 X 1023i^^,t/ V/^m-3 (H) 
^ ll.9xlO-ir.-/Xf cm- (He). 

Photons from the shocked stellar plasma diffuse into 
the jet through Thomson scatterings. Their penetration 
depth (mean-free path) in the shocked plasma frame is 
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The density of diffused photons inside the jet head drops 
off as 

n^{x) « n^^hS-''^^^'' (6) 
where x is the distance from the jet head radius Tj^t 0- 

III. PARTICLE ACCELERATION 

We calculate proton acceleration at the internal shocks 
which take place at a radius smaller than the termina- 
tion shock radius {rg^ < Tj^t)- In principle protons can 
be accelerated at the termination (forward and reverse) 
shock too, however, inverse Compton losses significantly 
reduce the maximum proton energy there . 

The maximum energy up to which protons are accel- 
erated in the internal shocks is limited by synchrotron 
losses in the jet magnetic field. (Inverse Compton losses 
are in the Klein-Nishina limit, [Eq. H17|) ]. and less impor- 
tant than synchrotron). We estimate the magnetic field 
Bjet in the internal shocks, in the jet frame, from the 
relation ^^:a^rl^cT'^^^^B'^^^JS■n = ^bL\^°- Dropping a - 1 
for notational simplicity and assuming a magnetic field 
equipartition fraction — lO^^^^-i, we get 
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The number density of protons in the stellar plasma 



frame is Ur, = 2r, 
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Equating the synchrotron cooling time 
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to the acceleration time 
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(both evaluated in the jet frame) we get the maximum 
proton energy in the observer's frame as 



30/100 ^ 0.3 s, for bulk Lorentz factor F = 100, is much 
longer than the acceleration time in Eq. ^ for all but the 
highest proton energies achievable. Hence the jet cross- 
ing time does not introduce any additional constraints 
for protons to accelerate to the maximum observer-frame 
energies of Eq. H1U|) . leading to a similar upper limit. 



IV. PROTON INTERACTIONS 
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This is smaller than maximum proton energy achievable 
in the observed optically thin GRBs, since the internal 
shocks take place at larger radii in the latter case. 

Electrons are also accelerated in the internal shocks, 
and the average electron energy in the jet frame is 

<7e>jet= (mp/TOe)£e ~ 3.7 X 10%. 2, whcrC Eg = 0.2£o.2 

is the fraction of thermal energy which goes into elec- 
trons. The corresponding peak synchrotron photon en- 
ergy in the jet frame is {3/ A){Bjet/ Bq) <je>%t "^e, 
where Bq = 4.4 x 10^"^ G. However, the jet is optically 
thick to Thomson scattering, with an optical depth in 
the jet frame, from Eq. of 



High energy protons accelerated in the sub-stellar in- 
ternal shocks undergo pj and pp interactions both inside 
the jet and near its termination boundary with the stellar 
envelope. We discuss each case below. 

Protons first interact inside the jet with synchrotron 
photons from co-accelerated electrons, and with the frac- 
tion of jet protons which are not accelerated at the inter- 
nal shocks. The dominant inelastic channel is P7 — > A 
with w 5 X 10~^^ cm^ cross-section (crp^_,A) and « 20% 
inelasticity. The P7 optical depth at the A resonance 
from Eq. (|13(l in the jet frame is then 
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Note that these values for sub-stellar jets are larger than 
those in the internal shocks responsible for 7-rays outside 
the star, due to the smaller distances and higher densi- 
ties. Hence the synchrotron photons thermalize to an 
approximate black-body temperature of 



■^7 jet 



can)] 



1/4 



n 9 r^/*r^^^^v 

U-^ ^52 ' 12.5 ^ 100 



1/2 



1.2L 



l/4^-l/2-p-l/2 



52 ' 11 



KeV 
^KcV 



(H) 
(He) 



in the jet frame. The corresponding number density of 
synchrotron photons in the jet frame is 



very large in both cases (see e.g. 7] for photomeson 
interactions in a different, optically thin, GRB acceler- 
ation scenario). With an optical depth so large, high 
energy protons lose all their energies through pry inter- 
actions. Secondary neutrons produced in the subprocess 
P7 — > A — > rnr^ also undergo photo-meson processes just 
like protons. This scenario is different from that of opti- 
cally thin shocks in observed GRBs, where the P7 optical 
depth is much smaller because of larger shock radii ]3,@1- 
The threshold proton energy for pr/ interaction at the A 

0.3r,et GeV2, 



resonance, from the condition E^'^E^^^^^ 
(12) in the observer's frame, from Eq. H12I) . is 
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We have also calculated energy losses by the protons 
due to photopair and photomeson interactions with syn- 
chrotron photons, which could possibly limit their accel- 
eration. The Bethe-Heitler process pj — > pe+e~ has a 
very large cross-section of « 5 x 10~^^ cm^ at threshold, 
but the energy lost by the proton is negligible (although 
they might play a role in more accurate calculations of 
the asymptotic energy losses). Photomeson interactions 
affect proton acceleration at energies only above A pro- 
duction threshold, which we discuss in the following sec- 
tion. The typical jet crossing time in the comoving frame: 



approximately an order of magnitude below the maxi- 
mum proton energy given in Eq. (|10|l . Due to the smaller 
(sub-stellar) radii of these shocks and the different target 
photons, the optical depths and threshold proton ener- 
gies are respectively higher and lower than those 4] in 
7-ray producing shocks outside the stars. 

Next, we consider pp interactions in the internal 
shocks. Typically a fraction < 1 of the protons are 
Fermi accelerated in the internal shocks of the GRB jet. 
The remaining fraction (1 — Cp) of cold protons pro- 
vide targets for pp interactions. The mean total cross- 
section for pp interactions in the TeV-PeV energy range 
is <crpp>« 6 X 10~^^ cm^. The corresponding optical 
depth in the jet frame, from Eq. is 

^'^pp^jet — (1 Cp^^pjet ^^pp^ fsh/^jet 
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much smaller than p7 optical depth. At energies well be- 
low the maximum proton energies, the competing effect 
of inverse Compton (IC) scattering reduces the pp scat- 
tering rate significantly, because of the large synchrotron 
photon density in the internal shocks. 

The cross-section for inverse Compton is ctic ~ 
aTh{me/mp)^ in the non-relativistic Thomson limit. The 
Thomson limit applies when E-yEp/nip <C 1 in the co- 
moving frame. Solving for the synchrotron photon ener- 
gies [Eq. ((T^ ] in the ri2.5 and ru cases, the Thomson 
limits are valid for El^^^^ < 4 x 10^ GeV and 7 x 10^ GeV 
respectively in the jet frame. The same quantity in the 
observer's frame 
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energy E^^^ Eq. ^ to fall below E^J^^J^. Solving 
Eq. 1(201) for X and t^^{x) = 5, we have x/l^ « 8. The 
shocked photon density at this distance, from Eq. (0J, is 
then 



< « 1.4 X lO-i^r-f r-/^m-3 



(H) (21) 



in the jet head. 

Protons below E^\^ in the ri2.5 case, interact with cold 
stellar protons in the H envelope. The number density of 
stellar protons is rip.H ~ 6 x 10^^ cm^'^. The correspond- 
ing pp optical depth for cold stellar protons in the stellar 
plasma frame is 

<Tpp>*=<crpp> np,H(r* - r^et) w 2.5 x 10* (22) 

which is very large and forces all protons to undergo pp 
interactions. 



is higher than E'^l^ [Eq. The ratio of IC to pp 

optical depth for the two internal shock radii is 
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The proton flux is suppressed according to Eq. ifT^ be- 
low E}p for pp interactions in the internal shocks. In the 
ultra-relativistic Klein-Nishina limit: E^Ep/m^ ^ 1, the 
pp optical depth becomes larger than the IC optical depth 
with increasing proton energy. However p^ interactions 
are dominant above £'p^f^ anyway. 

Protons below E^^^, in the ru case, undergo pp in- 
teractions in the internal shocks with reduced flux [Eq. 
(|TS|l ]. In the ri2.5 case, most protons below E^^^^ escape 
the internal shocks to interact with shocked photons in 
the jet head and cold stellar protons. We discuss these 
interactions below. 

The threshold proton energy at the A resonance for pj 
interactions with shocked photons in the jet head, from 
Eq. ®, is 
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in the observer's frame. The corresponding optical depth, 
in the jet head frame, depends on the photon number 
density given in Eq. © as 
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Since protons lose 20% of their energy in each p^ inter- 
action, it takes about 5 optical depths for a proton of 



V. NEUTRINO PRODUCTION 

High energy neutrinos are produced both in pj and 
pp interactions dominantly through secondary pion (tt^) 
decays as tt^ jiiy^ — > ev^Vpi^fj,. Each of the secondary 
leptons share roughly 1 /4 of the initial pion energy. 

Because of the high p7 optical depth in the internal 
shocks [Eq. H14I) ]. pions and muons are produced in the 
thermalized synchrotron photon bath. The high mag- 
netic field and synchrotron photon density in the shocks 
may force high energy pions and muons, from p^ inter- 
actions, to lose their energies through synchrotron radi- 
ation and inverse Compton scatterings before they de- 
cay to neutrinos. The same is true for p7 interactions 
with photons from the shocked plasma in the jet head 
for the ri2.5 case. Neutrinos from pp interactions tak- 
ing place outside the jet head in the unshocked stellar 
plasma, however, are not affected because of the much 
smaller magnetic field and photon density in the plasma. 
We discuss synchrotron and inverse Compton losses in 
the GRB jet below. 

The synchrotron loss time for pions and muons in the 
jet frame can be calculated using formulae similar to 
Eq. © for protons. Equating this time to the parti- 
cle decay time in the jet frame, T^jet = 2.6 x 10~*77r,jet 
s and T^jet = 2.2 x 10~^7^jet s, we get the maxi- 
mum synchrotron break energies, in the observer's frame: 

psb _ max ac 
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The inverse Compton (IC) losses are more severe for 
muons than for pions. The IC cooling time for a parti- 
cle of mass rrii in the Thomson and Klein-Nishina (KN) 
limits is 



VI. NEUTRINO FLUX CALCULATION 

The proton energy distribution in the sub-stellar inter- 
nal shocks, in the observer frame, is given by 
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respectively in the jet frame. The ratio of IC cooling time 
to the particle decay time, for synchrotron photons from 
Eqs. l^&EIl: in the Klein-Nishina limit, is 
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where C,p is the injection fraction of protons into the ac- 
celeration process and Eop is the fraction of jet energy 
dissipated in sub-surface shocks (whose order of magni- 
tude is ~ 1/2). The neutrino flux, the same for v^^ 
and in the pp and P7 interactions, from a single GRB 
buried jet at a distance D is flO] 
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for pions and 
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for muons. As a result, above 'fTr-^^jet ~ mTr-fj./Ejjet, 
only « 26% of the pions produced by P7 interactions 
in the internal shocks, in the r'12.5 case, decay to i^^. 
Production of Ug and D^, from /i-decay, in the ri2.5 
case is suppressed by IC losses. All flavors of neu- 
trino production in the internal shocks from pj inter- 
actions is suppressed in the rn case by IC losses. Below 
77r;/^jet ~ 'Ti7r;/i/-E'7 jet neutriuo production is suppressed 
as E~^ due to IC losses in the Thomson limit. Simi- 
larly, for pp interactions in the internal shocks, neutrino 
production is also suppressed by IC losses. 

Neutrino production with the shocked photons in the 
jet head, in the ri2.5 case, is also affected by inverse 
Compton losses. The decay length cTtt;^ in the shocked 
plasma frame exceeds the shocked photon penetration 
depth [Eq. ®] for 7^ > 9 x 10^ and 7^ > 10^. As a 
result pions and muons produced in the pr/ interactions 
"see" an increasing photon density according to Eq. ((HJ . 
For simplicity we compute the ratios of IC loss times to 
the decay times for pions and muons, in the Klein-Nishina 
limit: 7,7 ^ 9 x lO'* and 7/i > 6 x 10"*, as 
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using Eqs. As a result, all flavors of neutrino 

production is suppressed by IC losses in the jet head. 
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from 



where fpp ~ min(l, <rpp>) and = min(l,Tp-y) 
Eqs. ijini&IIII) respectively. We have used the PYTHIA 
6.2 event generator widely used in high energy parti- 
cle physics, to simulate pp interactions. For our problem, 
the neutrino multiplicity, through pion decay, in the pp 
interactions can be written as 
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for each flavor of neutrino. However, observed neutrino 
fluxes, for different flavors, are greatly affected by the 
environments at which they are produced. We discuss 
our flux calculation method next. 

In the rii (He) case, all flavors of neutrinos produced 
from p7 and pp interactions are heavily suppressed by 
inverse Compton losses both in the Thomson and Klein- 
Nishina limits [Eqs. (HH&EH]. We have calculated 
the muon neutrino flux (see Fig. ^) from pp interactions 
below Ep^^ [Eq. H15|l ]. in the internal shocks, using ~ 
40% of the proton flux [Eq. ^] according to Eq. (jTSJ. 
Above Ep^^ (corresponding to E^ ^ 10^-^ GeV) we used 
the full proton flux to calculate neutrino flux from p7 
interactions (not shown in the figure). 

The ri2.5 (H) case is more complicated. 
(Eq. 221, corresponding to E^, ~ 10"' '^ GeV), we cal- 
culate neutrino flux through pp interactions with cold 
stellar protons in the H envelope. We used ^ 70% of 
the proton flux [Eq. (|29|l ] escaping the internal shocks 
according to Eq. H16|) for this calculation. Note that 
neutrinos thus produced in the H envelope are not sup- 
pressed by inverse Compton losses. Above E^'l^^ and be- 



Below E;% 



Eq. p5|l ] Vp, production from p7 interactions, 



low E^l 

at the jet head, is suppressed by IC losses [Eq. (^5)) ]. 
Similar IC suppression [Eqs. (|26l & I27|) ] applies for z/'s 
from p7, in the internal shocks, above E^^^^. An extra 
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reduction of the flux happens above [Eq. [221 due to 
synchrotron radiation by pions. 
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EJIoglO(GeV)] 

FIG. 1: Diffuse muon-neutrino flux E^^uEap , shown as solid 
lines, from sub-stellar jet shocks in two GRB progenitor mod- 
els, H (ri2.5) and He (rn), each for two sets of shock/jet radii 
and overlying envelope masses (with similar curves for v^, 
and Vt). These neutrinos arrive as precursors (10-100 s be- 
fore) of 7-ray bright (electromagnetically detectable) bursts, 
or as the sole signal in the case of 7-ray dark bursts. Also 
shown is the diffuse neutrino flux arriving simultaneously with 
the 7-rays from shocks outside the stellar surface in observed 
GRB (dark short-dashed curve); the Waxman-Bahcall (WB) 
diffuse cosmic ray bound (light long-dashed curves); and the 
atmospheric neutrino flux (light short-dashed curves). For a 
hypothetical 100:1 ratio of 7-ray dark (in which the jets do not 
emerge) to 7-ray bright collapses, the neutrino fluxes would 
be 100 times higher than those plotted here. 

The diffuse i'^ flux levels are shown in Fig. as solid 
curves for the H (upper two curves) and He (lower two 
curves) star models. The nominal set of jet and shock 
radii discussed in and above (a ^ 1, H: r\i.^ and 
He: rii) give the uppermost solid curve of each set. For 
comparison, we also show the flux for a different choice of 
jet and shock radii, plotted as the lower of each set of solid 
curves; these latter have rjet,H ^ 10^^ '^ cm, r^/i^ ~ 10^^ 
cm in the H case, and rjet.Hc 10^"'^ cm and r^/i Hc ~ 
lO^*''^ cm in the He case respectively. Note that, because 
of a component, the neutrino flux we calculated is 
higher in the < 100 TeV range compared to optically 
thin shocks As a result, one gets more events in 

a neutrino detector than one would from optically thin 
shocks [T^. 

The effects of absorption by vN interactions in the 
overlying stellar envelope at the highest energies plot- 
ted in Fig. ^are calculated using the cross sections 
and the overlying grammages above the jet head. The 
absorption is negligible in the general case (a ^ 1, up- 
permost solid curve in each model), which have roughly 
~ O.IM0 (over 47r) of overlying envelope. For jet and 
shock radii occuring further in, as in the lower of each set 
of curves, the overlying envelope material that we took 
is more massive, \Mq. In this case the absorption 
effects can become noticeable, especially in He models. 



and this might serve as a diagnostic of the sub-surface jet 
depth (or the overlying material) . For a solar mass above 
the jet head (distributed over 47r), the optical depth for 
vN interactions becomes larger than unity at energies 
> 2.5 X 10" GeV in the H case and E^ > 2.5 x 10^ 
GeV in the He case (seen as a cut-off in the lower set of 
He curves). Typically 80% of the i/-energy is transferred 
to the secondary lepton: e or /Lt. These /i's further decay 
to J^'s, however, they are subject to severe synchrotron 
and inverse Compton losses and do not contribute sub- 
stantially to lower energy flux. 

The diffuse fluxes are based on a conservative estimate 
of the source density from the observation that ~ 10'^ 
GRB are known to occur per year over the entire sky, 
based on 7-ray detections. If these are attributed to mas- 
sive stellar collapses, each GRB should be preceded by 
neutrino percursor events, starting tens of seconds (typ- 
ically ~ 30 s Q) before the 7-rays. The total j/^ diffuse 
flux from all of these (7-ray quasi-coincident) events is 
given by the solid curves in Fig. ^ (Considering mass- 
mixing along the way to Earth from 2; ~ 1, the ratio 
of flavors as observed from earth should become unity). 
The number of neutrino bursts and correspondingly the 
diffuse neutrino flux may be up to a factor lOOg higher 
than what is shown in Fig. ^ if the ratio of the number 
of buried jets which do not emerge (7-ray dark bursts) 
to that of those which do emerge is lOOg, where <? ^ 1. 

The number of muon neutrino events, in a km^ de- 
tector (IceCube (l^ e.g.), from a single GRB at redshift 
z — I with buried jet of duration At « 30 s is 0.6 and 
0.003 in the ri2.5 (H) and in the rn (He) cases, respec- 
tively, in the 1-100 TeV energy range. The number of 
events can be signiflcantly larger for brighter {Ll^° > 10^^ 
ergs/s) or nearer (z < 1) bursts, or also if they are of 
longer duration {At > 30 s). 



VII. IMPLICATIONS 

The main feature of the high energy neutrino spec- 
tra of massive stellar collapse buried jet models related 
to GRB (and also of models involving a precursor su- 
pernova, Refs. 0, 01), is that a "thicker target" of 
nucleons and photons is available to give higher fluxes 
of ~ 1-100 TeV neutrinos, compared to those expected 
from internal or external shocks which occur outside the 
stellar progenitor 

The neutrino flux expected in association with electro- 
magnetically detected GRB, shown in Fig. ^ precedes 
the 7-rays by about 10-100 seconds, approximately the 
time taken by the jet to emerge from the collapsing core 
and the stellar envelope. Neutrinos coming from both 
P7 and pp interactions in the sub-stellar jet arc heavily 
suppressed by inverse Compton scatterings due to large 
photon density. The precursor neutrinos in the 1-100 
TeV range from a single 7-ray bright event, coming from 
pp interactions in the H envelope (ri2.5 case) can be de- 
tected with ICECUBE. The number of events from a sin- 
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gle burst at a distance D = lO^*'^ cm (z ^ 1) and dura- 
tion At = 100 s would be ~ 4 in the ri2.5 H-envelope star 
case, compared to < 0.04 in the rn He star case. (Rare 
single bursts occuring about once per year at D ~ 10^^'^ 
cm could yield ~ 10^ higher event rates per burst). Thus, 
detection of many bursts with 1-100 TeV neutrino pre- 
cursors in correlation with a subsequent 7-ray detectable 
GRB would be a strong indication for an H progenitor. 
A rare detection of a weaker 100 TeV neutrino precursor 
would favor a He star progenitor. 

The number of stars more massive than about 30 so- 
lar masses out to z ~ 1 undergoing a core collapse that 
leads to a black hole is about 10^/year corresponding to 
^ 3 X 10^^ yr~^ galaxy"^ of type II supernovae rate. An 
interesting possibility ^ is that many or perhaps all of 
these lead to jets as in GRB, in which case most of these 
jets do not break through the stellar surface, and hence 
are not observed in 7-rays (at most ^ 10^/year can be 
observed as GRB). If the mean jet solid angle subtends 
1% of the sky, this means < 10^/year such 7-ray dark 
(choked) bursts whose jets are directed at earth. These 
would be detectable only through their neutrino emis- 
sion, which should be the same per burst as for the 7-ray 
bright bursts, but they would be up to 100 times more 
numerous. The diffuse i/^ flux from these 7-ray dark col- 
lapses would be similar to that of the 7-ray bright solid 
neutrino curves of Fig. ^ but up to 100 times higher. 

The diffuse neutrino signals from both 7-ray bright and 
dark collapses should be detectable with IceCube, as can 
be seen from the sensitivity curves in Fig. ^ For the 
7-ray bright case plotted, the TeV-PeV flux detection is 
aided by the positional and temporal coincidences with 
electromagnetic detections. The typical angular resolu- 
tion of planned neutrino telescopes at TeV energies Q 
should he 6 < 1°. For the 7-ray dark cases, the dif- 
fuse flux for a total number of collapses < 10^/year out 
to z 1 would be sufficiently above the atmospheric 



neutrino background, especially in the PeV range, that 
detection appears possible even without coincident pho- 
ton flashes. (The 7-ray dark collapses could, however, be 
possibly associated with faint supernova- like optical/IR 
events, weeks to months later, if the stellar envelope is 
ejected). 

The non-detection of a diffuse signal would be signif- 
icant only if we have direct evidence for neutrino pro- 
duction in GRBs. This may come from a detection of 
either the GRB/afterglow neutrinos, or from a detection 
of the stellar pp neutrinos in the H case. In this case, 
an upper limit on the diffuse flux can serve to constrain 
the rate of "dark" bursts and the progenitor nature. If 
there is no direct evidence for neutrino production, the 
non-detection of a diffuse signal may imply no protons in 
the GRB jet and/or no proton acceleration. 

The flux may exceed the WB bound 9] if opaque colli- 
sions occur where the pp {p^) optical depth is very high, 
which implies that the nucleons do not escape. In this 
case, the proton energy generation rate of these "hid- 
den sources" may exceed the energy generation rate in- 
ferred from cosmic-ray observations, and hence the neu- 
trino flux may exceed the WB bound. 

The neutrino flux level and spectra in detected GRBs 
would serve as a diagnostic of the outer dimensions of 
the stellar progenitor at the time the explosion occured, 
providing useful constraints on the identity and evolu- 
tionary scenarios of the progenitors, as well as provid- 
ing clues and a consistency check for possible association 
with supernova-like objects. The measurement of such 
neutrino bursts from 7-ray dark collapses would provide 
useful constraints on the total rate of massive stellar col- 
lapses at high redshifts, a quantity of significant interest 
for cosmological reionization, structure formation and in- 
tergalactic metal enrichment scenarios. 
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